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Abstract 

The Fourier spectral analysis of this object suggests a nonlinear pulsation of period w 14d 
with a 3:2 frequency lock between two vibrational modes. We discuss an extensive search that we 
have performed in the L, M, T^jf space for models of period 14-16 days that exhibit this particular 
resonance. We find models satisfying the resonance constraints both in the Cepheid and in the post 
AGB regimes. The nonlinear pulsations of several candidates are presented. The nature of the star 
remains a mystery as none of the models satisfies all the observational constraints. 



1 Linear models 

We have performed an extensive search of linear nonadiabatic radial models with periods of P=14-16 
days. We have chosen this value because nonlinear and nonadiabatic effects tend to reduce the period 
somewhat, but the results are insensitive to the exact value of the period. 

The hydrostatic models and their linear stability analyses are computed for composition parameters 
A=0.35-0.90, each with Z=0.004, Z=0.01 and Z=0.02. The OPAL @ and Alexander & Ferguson [| 
low temperature opacities are used. Convection is ignored. For a fixed value of the temperature and 
the luminosity we iterate the mass to get the desired pulsation period of ~ 15 days for the fundamental 
mode or for the first overtone. Repeating this procedure for different values of T and L we obtain the 
period ratio {Pi/Pq or P2/P1) as a function of these two parameters. 

We present the results for A=0.7 and Z=0.01 in Figure 1. Each curve has a fixed luminosity 
(L=2,OOOL0 for the most left and L=12,000Lq for the most right curve). The mass decreases from 
the left to the right in each sequence. The highest mass is between 3 and 6 Mq, on the right side, 
the calculations on the blue side of the sequences are terminated for masses 0.3 and 0.6 Mq when 
we encounter numerical problems due to the high luminosity to mass ratio. The dots represent the 
unstable (fundamental mode) models. Two groups of unstable stars are found. The first group 
around Tejj=5,000K is inside the Cepheid instability strip while the second one gives the region for 
the unstable post- AGB stars (for a general linear survey of PAGB stars see eg IH). 



p /p 

0.070 



0.065 
0.060 
0.055 



4000 6000 8000 10000 T 

eff 

FIG. 1. First overtone to fundamental mode period ratio as a function of temperature, luminosity ranging from 2000 to 
110001/0; dots denote linearly unstable models; on the left is the Cepheid group, on the right the PAGB group. 

Models with a 3:2 period ratio exist only among the Cepheids (L=2, 000-4,000^0) for this compo- 
sition of X=0.70, Z=0.01. However the observations indicate a higher temperature and a higher 
luminosity for this star than our search for these model parameters yield. 

Interestingly, the period ratio increases close to the 3:2 resonance, right were the instability occurs 
in the second, PAGB group. The ratio Pi/Pq reaches the highest value around rejj=9,000 K which, 
however, is still lower than the spectroscopic estimate. 
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FIG. 2. Period ratio (Pi/Po) as a function of temperature for different values of the hydrogen content {X). 
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FIG. 3. Period ratio (P2/-P1) as a function of temperature. L = 7000, 9000, 10000 and 13000 Lq (left to right). 



The dependence of the period ratio on the composition is illustrated in Figure 2. Here we present 
only the high temperature models for L=9,000Lq. The curves are given for the composition of X=0.9 
(top), 0.7, 0.6, 0.5 and 0.35 (bottom) with Z=0.01. There is a clear increase of the maximum value 
of the period ratio as the hydrogen content is enhanced. An increase to X=0.9 gives the right period 
ratio, but this is an unacceptably large X, especially since we are in the PAGB region. For X=0.5 the 
period ratio is still increasing at the blue end of the curve but both the mass and the linear growth 
rate become small. 

We have also performed tests with different values of metallicity. Different values of Z only bring 
about a small shift of the curves. We cannot push up the period ratio to reach the resonance by 
varying the metallicity in the case of normal X values. 

We now turn to first overtone pulsators. Here, there now also are solutions for the P2/Pi=2/3 
resonance for lower values of X. The period ratio as the function of the effective temperature is 
displayed on Figure 3. The curves are given for L=7000, 9000, 11000 and 13OOOL0 and for X=0.5 
(dotted lines) and X=0.7 (solid lines). The resonance occurs in the range Tejj=8000-9000 K. The loci 
of the resonance are slightly shifted to the higher temperatures for the X =0.5 case compared to the 
higher X models. 



2 Nonlinear models 



From the Cepheid instability strip we have selected a resonant model with the following parameters: 
M=3.8Mq, L=3081Lq, Teff=5156K, X=0.7 and Z=0.01 with Po=14.1d. The numerical hydrody- 
namics shows that the resonance does indeed lead to a phase-lock with alternating cycles (cf Q for an 
explanation of this phenomenon). The theoretical velocity and light curve are displayed on Figure 4. 
The light variation shows similar alternations as the EROS star, and the Fourier amplitude spectrum 
is very similar to the observed one. The amplitude of the variation (~ 1.8 mag) is however much higher 
than the observed value (~ 0.5 mag). These differences are perhaps not astonishing since both the 
temperature and luminosity of the model are much lower than the observationally estimated values of 
the star. 




We have also calculated nonlinear models of the PAGB stars close to the linear resonance. In 
general we find that the computed amplitudes of these models are very small (usually between 0.001- 
0.010 mag.) in agreement with calculations of PAGB models but in disagreement with the observed 
value for this star. In contrast to the Cepheid model the growth rates are very large and the models 
tend to undergo irregular (chaotic) pulsations. 



One of the few models for which the hnear frequencies satisfy closely the resonance criterion has 
the following parameters: M=I.U06Mq, L=9566Lq, T^ff =9300K, X=0.9 and Z=OM. For this 
model both the fundamental mode and the first overtone are linearly unstable (with the growth rates 
of 7/0=0.36 and 771=0. 37). However, for none of the PAGE models have we found any alternations 
similar to the observations. 

As far as overtone pulsators are concerned the linear results have indicated that the P2 to Pi 
resonance provides good candidates for modelling this unique star. However, the fundamental mode is 
also unstable, and all of our nonlinear models (initialized with the first overtone velocity eigenvector) 
converge to the fundamental mode pulsation after a short transient. 

3 Conclusion 

The nature of the P=14d phase locked star remains a mystery in spite of a thorough search with 

radiative linear and nonlinear stellar models. The linear study shows that both Ccpheid models 
and PAGE models have the right period ratio between the fundamental mode and first overtone. 
However, all of the Cepheid models are underluminous and too cold compared to the observations. 
Furthermore, numerical hydro computations show that the pulsation amplitudes for these models are 
too high. For PAGE models the resonance does not occur with reasonable composition parameters 
for the fundamental mode, but can occur for first overtone pulsators. However, our hydrodynamics 
models are not able to reproduce the observed alternations in the pulsations. The theoretical light 
curves of the high luminosity and high temperature models also have very low amplitudes contrary to 
the observations. 

It is important to get fiirther observational confirmation and a better estimation of the parameters 
of this unique variable star, because the modelling can be an important test for our understanding of 
the physics inside this kind of stars. 
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